Tension and solute depletion in multilamellar vesicles 
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We show that a metastable multilamellar vesicle ('onion'), in contact with excess solvent, can 
spontaneously deplete solute molecules from its interior through an unusual, entropy-driven mech- 
anism. Fluctuation entropy is gained as the uneven partition of solute molecules helps the onion 
relieve tension in its lamellae. This mechanism accounts for recent experiments on the interaction 
between uncharged phospholipid onions and dissolved sugars. 

PACS numbers: 82.70.Uv, 87.16.Dg, 61.30.St, 61.30.Pq 



The richness of self- assembled structures exhibited by 
amphiphilic molecules (surfactants) in solution has at- 
tracted intensive research Q. More recently, attention 
has been turned toward the behavior of amphiphilic 
structures under stress S. An example for such struc- 
tures are multilamellar vesicles ('onions') that form when 
a phase of surfactant bilayers is subjected to shear stress 
||. Onion structures pose intriguing questions related 
to metastability and flow H [|, ^ and have possible ap- 
plications, e.g., for encapsulation of solid particles and 
drug delivery They also appear in biology, e.g., as a 
key ingredient in the maintenance of the lung surfactant 
monolayer 0. 

Some properties of onions are rather well understood, 
e.g., the viscoelasticity of the onion phase S, dependence 
of size on shear rate 0, and kinetics of swelling upon 
shear-rate change [|| or dilution ||, |i| . Other important 
features remain unresolved, in particular, the mechanism 
of the dynamic lamellar-to-onion transition ]To| | . Here we 
address another aspect of onions that has emerged in re- 
cent experiments [0 — their curious behavior in the 
presence of dissolved molecules that do not directly as- 
sociate with the bilayers. We show that, as in the case 
of onion dilution [j5| , a key role is played by surface ten- 
sion appearing in the lamellae when the onion is under 
osmotic stress due to contact with excess solvent. The 
unique way in which tension affects the onion entropy 
may lead to a situation where molecules are depleted 
from the onion interior for purely entropic reasons. 

Typical onions are of micron size, whereas the inter- 
membrane spacing is usually of order 10 nm, and the 
thickness of an individual membrane is a few nm. Thus, 
a typical onion comprises a spherical stack of hundreds to 
thousands spaciously packed membranes. Since synthetic 
bilayers are usually symmetric and thus do not possess 
a spontaneous curvature |Q , the onion structure is gen- 
erally not a global free energy minimum; the equilibrium 
state of the membranes, depending on water content, is 
a flat lamellar (L a ) phase, a disordered 'sponge' (L 3 ) 
phase, or a vesicle (L 2 ) phase [jl], Q. Yet, the path to 
this ultimate state requires onion disintegration or coa- 
lescence, involving breakage of hundreds of membranes 
per onion, and may take days B. Thus, an onion is in 



general far from equilibrium with its environment. In- 
ternally, however, the stack may relax on shorter time 
scales through the occasional formation of passage de- 
fects between membranes || Here we focus on this 
intermediate time regime, which may last hours to days, 
where the onion is in quasi-equilibrium but cannot an- 
nex or release membranes. We ask how such a structure 
affects the distribution of small solute molecules in its 
vicinity. 

In uncharged systems, or when electrostatic interac- 
tions are screened by salt, lamellar phases are stabi- 
lized by the Hclfrich interaction jlij - an entropic, 
fluctuation-induced repulsion between membranes. For 
tensionless membranes the free energy per unit area of 
this steric effect is f s t(d) = bT 2 /(nd 2 ), where d is the 
inter-membrane spacing, T the temperature in energy 
units, k the bending rigidity of the membranes, and b a 
numerical prefactor whose value is under debate; analytic 
calculations [11 4| give b = 3?r 2 /128 ~ 0.2, whereas com- 
puter simulations ]TH] find b ~ 0.06. The steric repulsion 
is supplemented by van der Waals attraction and shorter- 
range repulsion coming from hydration effects ]l6| . These 
interactions determine the equilibrium spacing d, as well 
as the transition point at which the stack 'melts'. 

When a lamellar domain coexists in equilibrium with 
excess solvent, the inter-membrane pressure equals the 
external osmotic pressure. This fact was used in ex- 
periments to measure inter- membrane interactions [ [l7| . 
In general, however, a spherical onion in contact with 
solvent is out of thermodynamic equilibrium and an os- 
motic stress may appear. Although water molecules can 
permeate rather freely through the membranes on the 
relevant time scales, further onion swelling is hindered 
by the lack of additional surfactant molecules and the 
high barriers to membrane breakage. Mechanical equi- 
librium requires that any pressure difference between the 
interior and exterior be balanced by surface tension in 
the lamellae. Since tension can be transmitted between 
membranes through occasional passages, a tension profile 
should in general be established throughout the onion. 
Such tension profiles have been shown recently to sta- 
bilize onions against 'melting' in an excessively diluted 
environment H. 
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The fact that tension has a significant effect on mem- 
brane interactions in onions is demonstrated by a simple 
argument. The parameter that determines the extent 
of this effect is the ratio between the inter-membrane 
spacing d and the capillary length I = (2T/7T0-) 1 / 2 , a 
being the surface tension pl|. If a; = d/l <C 1, ten- 
sion is negligible and the steric interaction is given by 
Helfrich's tensionless expression. When x > 1 mem- 
brane fluctuations are strongly suppressed, leading to 
an exponential decay of the inter-membrane repulsion 
p8[ . Suppose, by negation, that tension were negligi- 
ble, x <C 1. The interior pressure then would be roughly 
p ~ 2b'T 2 /(nd 3 ), where the numerical factor b' is a cer- 
tain fraction of b, i.e., of order 0.01-0.1. If there is no 
external pressure, the interior pressure is balanced solely 
by tension, a = \Rp — b'T 2 R/(Kd 3 ), where R is the 
onion radius. This gives x 2 ~ b'(T / n){R/d). Since typ- 
ically R/d ~ 10 2 -10 3 and K ~ 1-10T, one finds x > 1, 
contrary to the initial assumption |l9| . 

We begin the analysis by recalling the effect of ten- 
sion on the thermodynamics of a single membrane in 
a lamellar stack. The elastic Hamiltonian governing a 
single membrane is written in Fourier space as Ji = 

where /i q is a single mode 
The fluctuations can be di- 



i/d 2 ^ 4 + oq 2 )\h^ 
of height fluctuations. 



vided, using the Helfrich patch size £ 14 , into long- 



wavelength modes, q < 7r/£, that 'sense' the confinement 
due to the adjacent membranes, and short-wavelength 
ones, q > 7r/£, that are confinement-free. The former 
modes give rise to a steric repulsion for which, in the 
case of nonzero tension, there is no exact expression. 
A self-consistent approximation |2(| yields f s t(d,a) = 
(bT 2 /nd 2 )(x/ sinhx) 2 , which adequately coincides with 
the tensionless expression for small x — d/l and with 
renormalization-group calculations [fl8| for large x. Since 
tension suppresses the steric repulsion, this contribution 
to the free energy decreases with tension. Yet the dom- 
inant effect of tension is to increase the total free en- 
ergy by reducing the fluctuation entropy of individual 
membranes. Being independent of the spacing d, this 
contribution, denoted here is usually ignored. (It 
is implicitly included in the surfactant chemical poten- 
tial.) In the current study, however, we need to consider 
it explicitly. It is calculated by tracing e~ n / T over the 
short- wavelength modes |q| > 7r/£. The resulting inte- 
gral, f a = irT f*!^ gdgln(l + a / nq 2 ) , is dominated by the 
microscopic cutoff 5, which is of the order of the mem- 
brane thickness. The result is 

7T 3 T S 2 (7 

/ CT = ^[(l + y)ln(l + y)-ylny], y = (1) 

Let us proceed to the thermodynamics of the entire 
spherical stack. We employ a phenomenological, 'Flory- 
like' approach |2lJ, where the steric effect is taken into 
account accurately whereas the other interactions are in- 
corporated in a second- virial term. Using the expressions 



for / st and f a we obtain the free energy density of the 
onion as 
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In Eq. (§) <f> = S/(d + 5) = N5 3 /(4ttR 3 /3) is the sur- 
factant volume fraction (N being the total number of 
surfactant molecules in the onion), and \ characterizes 
the net attraction balancing the steric effect. We here- 
after consider, for mere simplicity, uniform concentra- 
tion and tension profiles within the onion p3). Thus, 
the Gibbs free energy of the entire onion is written as 



G(T,p ,N,a) = (47r/3)i? 3 [/ 



Pol 



/ LirR 2 a, where po 



is an external pressure. 

We now apply the quasi-equilibrium assumptions men- 
tioned above. The mechanical equilibrium condition is 
obtained from the requirement (dG/dR)T. Po ,N,<T = 0, 
leading to a Laplace equation, p — po — 2a/ R = 0, with 



bT 2 



1 



x cosh x 
sinh 3 x 



-xF 



(3) 



In addition, the assumption of internal equilibration im- 
plies a single, uniform chemical potential, 



v(F<j) = s 3 (df on /dF, T , Pota 
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[(1 — F) sinh x + 2xcosh: 
(1 - F) 3 sinh 3 x 



(4) 



which, as expected, is an increasing function of a. 

Unlike systems that can only shrink in response to ex- 
ternal pressure, onions have the option to release tension 
instead. For example, if N and R are constrained so 
that 4> is fixed, one can still increase the pressure, as seen 
from Eq. (|^), making the tension decrease in response. 
One way to introduce external pressure in practice is to 
add large molecules {e.g., polymers) that cannot be ac- 
commodated in the inter-membrane spacings and thus 
remain in the surrounding solution, exerting an osmotic 
pressure on the onion |l7], |23|. The current Letter, how- 
ever, concerns solute molecules that can dissolve in the 
inter-membrane water layers. In this case the molecules 
will partition between the interior and exterior, and the 
question is whether their distribution will be even or not. 
Systems exhibiting depletion usually involve steric con- 
straints (as for the polymers mentioned above) or repul- 
sive interactions. We consider a case where neither of 
these is significant; the molecules can comfortably dis- 
solve in the inter-membrane layers while being neither at- 
tracted to nor repelled from the membranes. Here we en- 
counter again a delicate consideration of relaxation times. 
We assume that permeation of the molecules through the 
membranes is kinetically hindered, such that during the 
intermediate regime of interest they can be viewed as im- 
pinging on the membranes and exerting pressure. This 
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holds, e.g., for sugar molecules, whose diffusion through 
an onion stack may take days JO). One might wonder 
then how the distribution of molecules can relax. This 
indeed poses a severe experimental difficulty; it has been 
overcome, nonetheless, by techniques of freeze drying and 
rehydration |23j ] , where a relaxed distribution equivalent 
to several weeks of molecular diffusion can be achieved. 

Let us introduce a dilute reservoir of solute molecules 
of volume fraction ipQ. We add to the free energy of the 
onion a contribution G s — (4ir/3)R 3 f s , with 



0.06 



f s = (T/a 3 )[ip\n(iP/ip )-(iP-ip )} 



(5) 



where ip = Ma 3 /(AttR 3 /3) is the volume fraction of the 
solute molecules inside the onion (M and a being their to- 
tal number and molecular size, respectively). Naturally, 
G s by itself is minimum for an even partition, ip — ipo. 
Taking the variation of Gtot = G + G s with respect to 
R we get, however, an expected contribution to the pres- 
sure, p s = (T/a 3 )(ipo — ip), which generates a coupling to 
the onion variables a and <p through the Laplace equa- 
tion. We are thus required to find the value of ip that min- 
imizes the total free energy Gtot > given that the modified 
Laplace equation, p(<p, a) —p s (ip) — 2a/ R = 0, is satisfied. 
(For simplicity we take po = 0, i.e., assume that p s is the 
sole source of osmotic pressure.) In order to examine 
changes in the system (e.g., tension release) as the solute 
content ip is progressively increased, an additional global 
constraint must be supplied. One simple limit is to as- 
sume that the surfactant volume fraction <p remains fixed. 
In this case the surfactant chemical potential fi is bound 
to decrease with increasing ip (since a becomes smaller 
while <p is constant). Another simple limit is to constrain 
fi. In that case <p increases (i.e., the onion shrinks) with 
increasing ip to compensate for the decrease in a. (In 
practice, neither of these two constraints is expected to 
be strictly accurate; the system may choose another path 
over the uUp,<r) surface.) 

Figure [l] shows the total free energy per unit volume, 
g = G to t/(47ri? 3 /3), as a function of ip for a given external 
solute content keeping \i fixed. The parameter values 
(cp for i/j = 0, k/T, S/R, ipo) have been chosen to match 
the experimental system of Ref. 
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The value of x is 
taken close to the binodal point of the stack, x = 0.9xbin, 
Xbin = 0/(1 — <P) 3 , mimicking a situation where the bare 
onion is not so far from equilibrium coexistence with the 
surrounding solvent. As seen in the figure, the favorable 
solute distribution is uneven, with a considerable deple- 
tion of about 20%. 

The dependence of the favorable ip and the corre- 
sponding tension release on ip is depicted in Fig. ^A, 
demonstrating how the effect is enhanced as more solute 
molecules are added. Finally, we have plotted in Fig. 
the dependence on the attraction parameter %, showing 
that depletion is in fact enhanced as one approaches the 
binodal. The reason is that the tension in the bare stack 
becomes lower and the onion is thus more susceptible 
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FIG. 1: Dependence of free energy on solute volume fraction. 
The minimum is obtained for ip < ipo, i.e., solute depletion is 
favorable. Parameter values match the conditions of Ref. JLl| : 
(p = 0.59 (in the absence of solute), k/T = 10, S/R = 1(H, 
a = 8, b = 37r 2 /128, ipo = 0.05, x = 7.57 (10% off the binodal; 
see text). The free energy per unit volume is given in units 
of 2bT 2 /(K8 3 ). 



to the solute perturbation. This finding implies that a 
strong depletion effect can be achieved by slight devia- 
tion of the bare stack from equilibrium coexistence with 
the excess solvent. Sufficiently close to the binodal, the 
bare tension is so low that the depletion annuls it, and 
further increase of x produces no effect. For the param- 
eter values of Fig. |^B this occurs extremely close to the 
binodal. 

Intriguing experimental results have been reported 
recently by Deme et al. on the interaction of un- 
charged phopholipid onions with dissolved hydrocarbon- 
ates (small sugars) [jll). Using small angle neutron scat- 
tering and exploiting the distinct contrast match points 
in the presence and absence of sugar, they determined 
that the sugar volume fraction inside the onions was tens 
percent lower than its external value. In small angle x-ray 
scattering the peak corresponding to the lamellar period- 
icity was observed to broaden upon adding sugar, indi- 
cating an increase in lamellar fluctuations. The authors 
attributed the enhanced fluctuations to strong reduction 
in bending rigidity due to the added sugar, yet mentioned 
that there was no known mechanism for such an effect. 
Indeed, it is hard to see how molecules that are well dis- 
solved in the water layers and do not associate with the 
membranes could affect their rigidity. The mechanism 
suggested here provides a straightforward explanation for 
the increased fluctuations, as well as a relation to the ob- 
served depletion — the enhanced floppiness of the mem- 
branes is caused by tension release rather than reduced 
rigidity. 
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FIG. 2: Dependence of depletion (solid curve, left ordinate) 
and tension (dashed curve, right ordinate) on the external 
solute content (A) and attraction parameter \ (B). The 
depletion effect is enhanced as solute molecules are added and 
as x approaches its binodal value Xbin- Parameter values are 
the same as in Fig. |l]. Tension is given in units of T/S 2 . 



There is an important aspect that is not captured by 
the current model — the onions were found in the ex- 
periment to swell with increasing sugar content. One 
possible, previously studied swelling mechanism invokes 
the sugar effect on the dielectric constant of the solvent, 
which modifies the van der Waals attraction between 
membranes |l7j. (This would translate in our notation 
to a ^-dependent %.) Note that the depletion mecha- 
nism presented here is quite general and should apply, in 
principle, to other solute molecules such as dissociated 
counterions, provided that their permeation through the 
membranes is not too fast. We have considered, however, 
neutral molecules whose direct interactions are weak; in 
the case of charged molecules strong electrostatic inter- 
actions will set in and require a more complicated model. 



We have treated the interplay between various ingredi- 
ents of this complex system on a basic, phenomenological 
level. In particular, we have assumed all quantities (sur- 
factant and solute volume fractions, membrane tension) 
to be uniform throughout the onion which, in fact, can- 
not be strictly correct The next natural step would 
be to consider nonuniform profiles and see, e.g., whether 
the solute forms a nontrivial profile inside the onion. It 
would be interesting to check also whether the unusual 
depletion mechanism suggested here is exploited in bio- 
logical systems, such as the lamellar bodies in the lung. 
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